Magnetic field dependent impact ionization in InSb 
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Carrier generation by impact ionization and subsequent recombination under the influence of 
magnetic field has been studied for InSb slab. A simple analytic expression for threshold electric 
field as a function of magnetic field is proposed. Impact ionization is suppressed by magnetic field. 
However, surface recombination is dependent on the polarity of magnetic field: strengthened in one 
direction and suppressed on the opposite direction. The former contributes quadratic increase to 
threshold electric field, and the latter gives additional linear dependence on magnetic field. Based 
on this study, electrical switching devices driven by magnetic field can be designed. 
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PACS numbers: 

There are many reports about the impact ionization 
in semiconductors. However, it is hard to find the one 
explaining the influence of magnetic field on the impact 
ionization. To the best of our knowledge it is only re- 
cently that a reasonable pilot model has been proposed 
for the magnetic field effect and compared with exper- 
imental results^. The previous model is restricted to a 
special case: electric transport is quasi-ballistic, and car- 
rier recombination is independent of magnetic field. We 
approach to this issue from more general background, in 
which an electron can experience many scatterings before 
reaching the impact ionization, and magnetic field affects 
carrier recombination process. 

Our model gives a result that magnetic field con- 
tributes to the carrier generation and recombination pro- 
cess: the field reduces the generation rate and increases 
the threshold voltage, and it also makes recombination 
rate sensitive to the polarity of magnetic field. After de- 
scribing some picture regarding our model qualitatively, 
quantitative treatment will be followed. 

When a high bias voltage is applied, electrons acceler- 
ate to a high speed. If the kinetic energy acquired from 
the electric field equals the ionization energy, impact ion- 
ization occurs. Upon impact with the lattice, the electron 
expends its kinetic energy on ionizing a valence electron 
(refer Fig. 1 (a)). This process produces electron- hole 
pairs and abruptly increases the electric current in the 
device. Impact ionization makes equal number of ex- 
cess electrons and holes. Because electron mobility is 
more than 100 times larger than that of hole in InSb, we 
consider only electronic conduction in impact ionization 
regime. Before reaching the ionization energy, the ener- 
getic electron can experience energy loss due to inelastic 
scatterings. To achieve impact ionization, the electron 
should accumulate kinetic energy despite the inelastic 
scatterings. Magnetic field affects this carrier genera- 
tion process. When magnetic field is applied, the Lorentz 
force deflects the electronic trajectory, and the net gain of 
kinetic energy for a given path length is reduced (Fig. 1 
(c)). To achieve the ionization energy a longer trajectory 
is required; however, this longer trajectory gives rise to 



the greater possibility of inelastic scattering. Thus, the 
deflection of the electronic trajectory caused by magnetic 
field leads to suppression of the impact ionization. To re- 
store impact ionization, a greater electric field is needed 
to increase the net energy gain between the scatterings. 
Consequently magnetic field suppresses the carrier gen- 
eration and increases the threshold electric field. 

Recombination is an elimination process of electron- 
hole pairs and generally follows carrier generation. We 
are interested in recombination at the two interface, SI 
and S2, as depicted in Fig. 1(d). SI has higher recombi- 
nation velocity than that of S2 and carrier electrons are 
readily rccombined near SI. Magnetic field produces the 
Lorentz force. Electrons accumulate or deplete near SI 
according to the polarity of magnetic field. When the po- 
larity of magnetic field is negative, negative z-direction in 
Fig. 1(d), the Lorentz force deflects electrons to SI and 
recombination process is facilitated, whereas positive po- 
larity makes electrons near SI depleted and results in 
relatively slow recombination. 

When the bias voltage exceeds threshold voltage, gen- 
eration rate is larger than recombination rate, and then 
number of electrons increases with time, which is known 
as avalanche state. Under steady-state, however, genera- 
tion process is balanced with recombination, i.e., genera- 
tion rate is same as recombination rate. In this work we 
consider the limit of steady-state, which is on the border 
of avalanche state. 

Now we are in position to treat the model quantita- 
tively. The model proposed in this work presents a sim- 
ple analytic expression for the threshold electric field as a 
function of magnetic field. For small band gap semicon- 
ductors such as InSb, the ionization energy is approxi- 
mately equal to the band gap energy e^. Energetic elec- 
trons undergo inelastic scatterings before their kinetic en- 
ergies reach e g . The dominant scattering process of InSb 
at room temperature is optical phonon scattering 3 -^, and 
the optical phonon energy fnu is known to be 23 meV— . 
Hence, in the present model each scattering makes an 
energy loss of hw . We will obtain a probability for an 
electron to acquire e g in spite of energy losses due to in- 
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FIG. 1: (Color online) Schematic diagrams illustrating transport mechanism in impact ionization process, (a) An 
electron accelerated by electric field achieves its kinetic energy equal to ionization energy and generates electron-hole 
pair. Before reaching ionization energy, the electron experiences inelastic scattering, (b) Inelastic scattering can be 
classified to two groups: small- and large- angle scattering, (c) Between the scatterings, the trajectory of the electron 
is deflected by magnetic field B, which suppresses acquisition of kinetic energy and therefore impact ionization, (d) 
Recombination is dependent on the polarity of magnetic field when recombination strength at SI and S2 are different. 
Axis and two interface, SI and S2, in the sample are depicted in the bottom of (a). 



elastic scatterings. Assuming this probability is propor- 
tional to carrier generation rate, generation rate will be 
expressed in terms of electric and magnetic field. Intro- 
ducing steady-state condition and recombination param- 
eters, a simple analytic relation between the threshold 
electric field and magnetic field will be proposed. 

Adopting Dumke's theory for InSb-, scatterings are 
classified into two groups according to scattered direc- 
tions: small- and large-angle scattering (Fig. 1(b)). For 
an electron incident along the direction of an electric 
force, the small-angle scattering produces a scattered an- 
gle less than n/2, and the electron is ready to be ac- 
celerated again by the electric field after the scattering. 
The large-angle scattering gives a scattered angle greater 
than 7r/2, which results in deceleration and lose a chance 
to acquire further kinetic energy from the electric field. 
Thus, large-angle scattering should be avoided to achieve 
impact ionization. 

Our model starts from this classification. The large- 
angle scattering probability in a time interval of dt can 
be given by (M/tl, where tx is the relaxation time of 
the large-angle scattering. Then, the probability P 
of an electron surviving large-angle scattering is P = 
exp(— J l/r^di) for a finite period of time. P can be 
expressed with an electric field. An electric field E sup- 



plies energy to an electron at a rate of eE(v x ), where e 
is the electron's unit charge, and (v x ) is an average ve- 
locity parallel to the electrostatic force in the time inter- 
val between the successive small-angle scatterings, i.e., 
(v x ) = (1/t s ) J q 3 v x dt, where t s is the relaxation time 
of the small-angle scattering. Energy loss by the small- 
angle scattering is given by fiw /r s . Thus, the net rate 
of energy gain is expressed by de/dt = eE(v x ) — hio /T s 
. Using dt = de/(eE(v x ) — Hui /t s ), P is rewritten as 

p= exp (r--- , * , V a) 

\Jhu> a t l eE(v x ) — hio /T s J 

A reasonable choice for the lower bound in the integrand 
is huj because optical phonon scattering is absent for an 
electronic energy less than Hu) . When an electron with 
an initial energy of huj moves to a final energy states 
of Eg, various paths are possible, and P is dependent 
on these paths. We only consider a path which gives 
a maximum value of P. The path having the shortest 
electronic path- length gives the maximum P. A diffusion 
effect in the energy space is not considered; it is a second 
order effect^ and ignored in this work. 

To elucidate (v x ) in Eq. (JTJ, more details for the small- 
angle scattering are needed. Velocity is generally depen- 
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FIG. 2: (Color online) (a) The parameters in Eq. ([3]) ver- 
sus energy. The curve of t s /t l indicates that small angle 
scattering is preferable to large an gle scatter ing. Voo is an 
asymptotic velocity defined by yje g /(2m c ). (b) M ver- 
sus mobility times threshold field /i Eq . The magnitude 
of M and \iEq are restricted: 0.068 < M < 0.24 and 
/j,Eq > 4.7 x 10 5 m/sec. Note that only intrinsic prop- 
erties of InSb (e g , m c and Huj ) are used to obtain (a) and 

(b). 



dent on kinetic energy. However, due to the band non- 
parabolicity 8 of InSb, the velocity of InSb is a slow func- 
tion of energy for energetic electrons (refer Fig. 2(a)). 
Hence, the variation of the magnitude of velocity be- 
tween the successive small-angle scatterings is considered 
to be negligibly small even though the corresponding en- 
ergy change is considerable. After a scattering event the 
electron may have various scattered directions, and these 
directions work again as incident directions for the next 
scattering. The shortest path is achieved when this scat- 
tered direction is in parallel with the electric force (along 
x-axis). Thus, the direction of electronic velocity just af- 
ter the small-angle scattering is considered to be in par- 
allel with the electric force. 

Magnetic field is involved with (v x ) in Eq. (p}. In the 
presence of a magnetic field B, classical trajectory of an 
electron is governed by the Lorentz force. For an elec- 
tron with velocity of v moving toward the electric-force- 
direction (refer Fig. 1(c)), the classical trajectory gives 
l/r s Jp 3 v x dt — vsin(Lu c T s )/uj c T Sl where w c = eB/m* and 
m* is electronic effective mass. Using the approximation 
sin(uj c T s ) = 1 — 1/6- (lo c t s ) 2 for low magnetic field, (v x ) 
in Eq. ([T]) can be replaced with v (1 — 1/6 • (lu c t s ) 2 ). The 
total scattering rate l/r s + 1/tl is nearly independent 
of energy^, and t s is favored over T£ in the overall en- 
ergy range (Fig. 2(a)). Thus, electronic mobility \i is 
determined by r s , which allows t s — fim c /e. m c is the 
effective mass at the conduction band edge. We consider 
E as a small deviation from the threshold electric field at 
zero magnetic field, Eq. The right term in Eq. ([T]) can be 
expanded to the first order of {E/Eq — 1) — 1/6 • (uj c t s ) 2 . 




E=E1 



I E=E2 



FIG. 3: (Color online) Normalized threshold electric field 
E/Eq — 1 as a function of mobility times magnetic field 
fi B for a variety of ratios of surface recombination ve- 
locity, Si/s2- Eq. is used on the assumption that 
ftE = 5 x 10 5 m/sec, S2 = 10 4 m/sec and /i p //i is effective 
mass ratio of electron and hole, 0.03. The yellow area rep- 
resents avalanche state for s\/ S2 — 3. The avalanche state 
is achieved for an given electric field E = El (E = E2) 
when [iB> —0.5 {\iB > 0.5). This characteristic provides 
abrupt change of current in current versus \iB plot. The 
inset shows an expected schematic trace for E — El and 
E2, which can be used as an electrical switch controlled by 
magnetic field. 



Following some algebra, it is readily shown that 

P/P - 1 = C G ((E/E - 1) - (l/6)M(/iP) 2 ), (2) 
where 

C G 
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and Po is amount of P at E = E and 5 = 0, i.e., 
exp (- f' 9 a- P d£ r ) . 

The parameters in Eq. §3§ can be calculated as fol- 
lows. An analytic form of t s /t^ is given by Log((l + 
k r ) 2 /{l + kf))/Log((l + k 2 )/(l - k r ) 2 ) 6 , where k r = 
k(e) — hjj /k(e) and k(e) is a wave vector of an elec- 
tron for a given energy e. The non-parabolic disper- 
sion is given by e(k) = —e g /2 + \Js 2 /^ + e g h 2 k 2 /(2m c ) 

from which the following expressions can be ob- 
tained: k(e) = (y / 2m c /h)y/e 2 /e g + e, v = (l/h)de/dk 
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= (l/V2m c )^/£ 2 £g + ee 2 g /(e + e g /2) and m c /m* = 

(m c /h 2 )d 2 e/d 2 k = (2e/e g + 1)~ 3 . In InSb e g is 0.17 eV 
and m c is 0.013 times the mass of a free electron. Using 
these relations, Co and M in Eq. ([3]) can be obtained. 

It is only a few electrons which survive large-angle scat- 
tering and obtain kinetic energies equal to e g and finally 
contribute to the impact ionization 10 . Because the im- 
pact ionization is governed by the energy-gain process 
specified by P, a generation rate per unit carrier due to 
the impact ionization is assumed to be proportional to P. 
Therefore, P/Po — 1 in Eq. © represents a normalized 
generation rate. 

Recombination relies on magnetic field when the re- 
combination strength at the two interfaces are different 
(see Fig. 1 (d)). In low magnetic field regime the first 
order approximation of recombination rate R(B) with re- 
spect to B gives R(B)/R — 1 = —Cr[j,B, where R is 
recombination rate at zero field 11 . When surface recom- 
bination is dominant over bulk recombination and sam- 
ple thickness is smaller than carrier diffusion length, the 
following simple expression can be obtained from Lile's 
results^: 

Cr =(---) VpEo- (4) 

S\ and S2 are surface recombination velocities at the two 
interfaces, SI and S2, respectively, and difference of them 
gives non-zero value of Cr. The steady-state condition 
asserts that the generation rate is equal to the recombina- 
tion rate, which leads to P/Po = R(B)/R and therefore 
P/P - 1 = -Cr fj, B. Then, the threshold field in Eq. © 
is expressed by 

E/E - 1 = \M(nB) 2 - %^B. (5) 



Plotting threshold field E according to magnetic field 
permits an overview of the present model. The plots in 
Fig. 3 are calculated ones using Eq. (j4j and Eq. ([5]). 
SI unit for mobility and magnetic field makes /x B di- 
mensionless. For identical surface recombination at the 
two interface, s\ = S2, the normalized threshold field in- 
creases quadratically with magnetic field. This increase is 
caused by the suppression of impact ionization by mag- 
netic field. Difference in the surface recombination ve- 
locities, s± > S2, adds the linear term and the thresh- 
old field becomes asymmetrically dependent on magnetic 
field. Note that the positive (negative) polarity of mag- 
netic field corresponds to the middle (bottom) diagram 
in Fig. 1 (d). Weak (strong) recombination results in 
small (large) threshold electric field. 



The curves in Fig. 3 represent boundaries between 
normal and avalanche states in the space of electric and 
magnetic field: the avalanche and normal states corre- 
spond to the upper and lower areas of the curve, respec- 
tively. The yellow area corresponds to avalanche state 
for s\/s2 = 3, for instance. By varying magnetic or elec- 
tric field, one of the two conducting state, normal and 
avalanche state, can be selected. An interesting applica- 
tion of this phenomenon is switching device. For a given 
electric field, electric current can be changed abruptly by 
varying magnetic field, which can be a good candidate 
of magnetic-field-driven electrical switching device (refer 
the inset of Fig. 3). 
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